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A recently proposed "DFT+dispersion" treatment (Rajchel et al., Phys. Rev. Lett., 2010, 
104, 163001) is described in detail and illustrated by more examples. The formalism derives the 
dispersion-free density functional theory (DFT) interaction energy and combines it with the dis- 
persion energy from separate DFT calculations. It consists in the self-consistent polarization of 
DFT monomers restrained by the exclusion principle via the Pauli blockade technique. Within the 
monomers a complete exchange-correlation potential should be used, but between them only the ex- 
act exchange operates. The applications to wide range of molecular complexes from rare-gas dimers 
to H-bonds to 7r-electron interactions show good agreement with benchmark values. 



I. INTRODUCTION 

Density Functional Theory (DFT)-based methods pro- 
vide the most important viable approach to large systems 
of nano- and biotechnological relevance. Their success for 
determination of structure, energetics and other static 
properties of dense matter is well known. 

However, treatment of weak non-covalent interactions 
by DFT remains plagued by spurious and erratic results. 
This is a consequence of the fact that stabilization in 
these complexes is determined by dispersion interaction, 
not accounted for in standard DFT functionals.^'^ Con- 
sequently, exchange-correlation potentials derived from 
local and semi-local models and combined to satisfy em- 
pirical data often feature artifacts when applied to sys- 
tems with large non-local correlation effects. 

Serious efforts have recently been invested to incor- 
porate the dispersion effect into the DFT framework, 
and the results are promising. ^^^'^ At the same time, 
in many practical applications remarkable progress has 
been achieved by using a posteriori dispersion correc- 
tions, model and/or semi-empirical, added on the top of 
regular DFT calculations (DFT-|-D)^ in the spirit of the 
classic SCF-|-dispersion model of Ahlrichs et al.^"* and Wu 
et al.6 

In order to make a DFT-|-dispersion strategy success- 
ful, one needs two ingredients: a dispersion-free DFT 
interaction energy^ and the model for the dispersion en- 
ergy. The first, the DFT dispersion component, may 
be obtained via the SAPT approach, but also from 
other approximate DFT treatments. ^'^'^^'^^ The second 
ingredient, the dispersion- free DFT interaction energy, 
is commonly obtained as the supermolecular DFT inter- 
action energy. However, in contrast to the SCF inter- 
action energy, the composition of the DFT interaction 
energy is neither understood nor controlled. In particu- 
lar, an exchange-correlation functional is always an ap- 
proximation to some extent delocalized, and never ex- 
actly dispersion-free, and the problem of double count- 



ing of obscure dispersion terms arises. In addition, the 
exchange-correlation potential compromises the demands 
of many users and various training sets, and thus is prone 
to erroneous behavior. 

A rigorous approach requires a DFT interaction energy 
that a priori neglects non-local long-range interaction en- 
ergy terms (dispersion) but allows for accurate mutual 
exchange and mutual polarization effects — an analogue 
of the SCF interaction energy at the DFT level of theory. 
Such a DFT interaction energy could be confidently and 
rigorously supplemented with a dispersion component 
obtained also at the DFT level of theory using SAPT^^-^^ 
or other formalisms. ^'^'^^ The goal of this work has been 
to define an accurate DFT-|-dispersion treatment which 
is based the derivation of "dispersion-free" interaction 
energy arising between DFT monomers to which a pos- 
teriori DFT dispersion energy is added. The former, the 
dispersion-free DFT interaction energy, leaves out any 
intermolecular correlation terms — thus avoiding artifac- 
tual exchange and doubly-counted dispersion terms. The 
latter may be obtained from e.g. TDDFT and/or SAPT 
calculations. The basic formulation of the method and 
selected results for model systems: rare-gas dimers and 
H-bonded systems, have already been presented by Ra- 
jchel et al.^. 

To this end a novel hybrid DFT approach for calcu- 
lations of van der Waals complexes has been proposed 
by Rajchel et al.^. It uses the formalism of Gutowski 
and Piela^° termed Pauli blockade-Hartree Fock (PB- 
HF) combined with the bifunctional formulation of DFT 
of Rajchel et al.^^. At the first stage, the subsystems' 
sets of orbitals are separated by symmetrical orthogonal- 
ization of the orbitals residing at different monomers. At 
the second stage, one iteratively evaluates the interaction 
energy between two DFT monomers described by Kohn- 
Sham determinants in a manner analogous to the HF 
method. That is, the monomers are mutually polarized 
until self-consistency in each other fields is reached un- 
der the constraint of the Pauli exclusion principle (Pauli 
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blockade technique). Various exchange-correlation DFT 
potentials may be used within the monomers (see below) , 
but the intermonomer exchange-correlation potential is 
reduced to only the exact exchange, and thus neglects the 
dispersion contribution. In the third stage, the disper- 
sion component is a posteriori added, from SAPT(DFT) 
calculations. In such a way, the erratic behavior of ap- 
proximate exchange functionals around the equilibrium 
separations is no longer the case, while the missing long- 
range attraction caused by dispersion is added without 
the problem of double counting. 

It is worthwhile to note that the above approach bears 
resemblance to the range-separation idea in DFT^^"^"* 
differing in that it is based on separation of monomers 
rather than ranges. 

In Sec. II the "dispersion- free" PB(DFT) is derived as 
a particular case of the PB treatment of DFT reported 
by us recently.^^ In Sec III we report numerical results 
of our method for representative van der Waals and H- 
bonded systems, and discuss the overall performance of 
the method. 



II. THEORY 

A. Pauli blockade method 

Supermolecular energy in terms of DFT can be defined 
as the difference between the total energies of the dimer 
AB and the individual monomers A and B, separated to 
infinity: 



DFT 



^int — ^AB 



pDFT 



pDFT 



(1) 



The interaction energy calculated in this way contains the 
correlation contribution, however the long-range correla- 
tion effects are not taken into account correctly within 
standard density functionals. 

It was demonstrated by Gutowski and Piela'^" that 
the Hartree-Fock (HF) supermolecular interaction en- 
ergy may be exactly recovered by solving the HF equa- 
tions for monomers in the presence of the external per- 
turbation consisting of the electrostatic potential and 
the non-local exchange potential generated by the sec- 
ond monomer. They have also proposed a convenient 
computational scheme in terms of mutually orthogonal- 
ized A and B monomers' occupied orbitals, {a^j^g^ and 

\ hk \ (the quantities expressed in the orthonormal- 
L J fees 

ized orbitals are henceforth marked with tilde). With 
such orbitals, satisfying 



0, 



(2) 



the supermolecular interaction energy in Hartree-Fock 
method can be written as 



(3) 



where 



AEa = i?r - El 



(4) 



is the difference between the final and isolated 
monomer A HF energy. 



Eaut — 
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is the electrostatic interaction, 

PA(r) =2^|a,(r) 



ieA 



is the monomer A density, 

<(r) = - E 



aeA 



|r-R„| 



(5) 



(6) 



(7) 



is the monomer A nuclear potential with a labeling the 
coordinates of monomer A's nuclei, each described by its 
position and charge Z^, 



Mn^-EE 
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is the intermonomer nuclear repulsion energy (constant 
for a fixed geometry), and finally 



E 



-2^ ^ {a,bk 

ieA keB 
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(9) 



is the exchange interaction. The monomer A orbitals are 
the eigenfunctions of the following modified Fock opera- 
tor: 



/a + if' + 



(10) 



where /a is the standard Fock operator build of {aOigA 
orbital set, and the two remaining terms are electro- 
static and the non-local exchange potentials generated 
by monomer B, respectively. Monomer B orbitals are 
obtained analogously. 

The procedure introduced by Gutowski and Piela^° can 
be generalized to the case where the two subsystems are 
described by KS orbitals. To this end we note that the 
total density of the system can be represented as the 
sum of monomer densities obtained from orthogonalized 
orbitals: 



Mb = Pab = Pa + Pb, 



(11) 



where the monomer A density is calculated as in (6) but 
with orbitals being the solutions of the following modified 
KS equation: 

(/F(r) -t- A^r(r) + «B'^'(r))a,(r) = e^,Mv) (12) 
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and satisfying (2). In Eq. (12), /^^ is the standard KS 
operator built of {aijieA orbitals, 



is the non-additivity of the monomer A exchange- 
correlation (xc) potential, and v^^^ (r) is the electrostatic 
potential of Eq. (10). With pAB decomposed in such 
way, the total KS interaction energy may be written as 
a functional of pA and pe densities: 

Efn^ [pa; Pb] = AEa [pa] + A^B [pb] + 
+ E^ist [pa; Pb] + AE^c [pa; Pb] ■ (14) 

The terms of (14) are as follows. Monomer A deforma- 
tion is 



A^A [pa] = Ea [pa] - Ea [pI] 



Finally, the exchange-correlation interaction is calculated 
in a supermolecular manner, 

(13) AE^,[pa;pb]=E'''[pa+pb]-E^'[pa]-E^'[pb]. (23) 

Inserting monomer densities calculated using Eq. (6) with 
orbitals satisfying (2) and (12) into (14) one restores the 
DFT supermolecular interaction energy (1). The expres- 
sion (14) is potentially exact, i.e. it yields the exact in- 
teraction energy provided that the exact xc potential is 
used. 

Technically, the orbitals satisfying (12) are found in 
a self-consistent iterative process and the orthogonality 
condition (2) is imposed by the brute-force incorporation 
of the penalty operator and successive Lowdin orthogo- 
nalization. Depicting iteration numbers in square brack- 
ets, the nth iterative step for monomer A reads 



(15) 



and Pa the density of the unperturbed A monomer, 
i.e. the density built of orbitals {a°}^g^ satisfying un- 
perturbed KS equations. 



/-^■'^(r)a^(r) = 4,a^(r). 



The total energy of monomer A is expressed via the stan- 
dard KS functional, 

Ea [pa] = [pa] + [pa] + J [pa] + 

+ E^- [pa] + Vr- (17) 

The functional (17) includes the non-iteracting kinetic 
energy: 



[PA] = -Y.i^^\^r]~a.), 

nuclear-electron attraction energy: 

^a"'[pa]= / <(r)pA(r)d^r, 

Jm.3 

coulombic energy: 
1 



(18) 



(19) 



J[pj, 



PA(ri)pA(r2)^3r^^3,^^ (20) 



2 Jr3 Jr3 ri2 
exchange-correlation (xc) energy: 

[-^a] = ^3 ^'''(pA(r); {VrPA(r); . . . }) d\ (21) 

which is evaluated through the numerical integration of 
the F^^ integrand on a grid of points around monomer A. 
The last term of (17) is the monomer A nuclear- nuclear 
repulsion energy. Similar expressions can be written for 
monomer B. The electrostatic part of the interaction en- 
ergy has the same form as in Eq. (5) and in terms of 
densities it can be easily rewritten as 



Ecist[pA; Pb] = 

= I vl%v)pA{v)d'r+ f vT{rypB{r)d' 
ri2 



r+ 



(22) 



' ?KS[ri-l] 



A^xcfn — 1] 
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(Ig) where the penalty operator is 



keB 



(24) 



(25) 



and it is obvious that its action on monomer A's occupied 
orbitals annihilates them once the orbitals are orthogo- 
nal. ?7 > is a scaling parameter not affecting the final 
solutions. The equivalent of (24) for monomer B is ob- 
tained through the interchange of the A and B subscripts. 
The orbitals obtained in each iteration are orthogonal- 
ized, yielding an orthonormal 



keB 



(26) 



set. The iterations start with the unperturbed orbitals 
obtained in Eq. (16) and its analogue for monomer B. 

The zeroth order interaction energy may be viewed 
as an analog of the well-known HF-based Heitler- 
London (HL) interaction energy. Specifically, we define 
the DFT-based HL interaction energy as 



e^ = Eab [pI;p"b]-e[pI]-e[p'^], 



(27) 



where the densities Pa Pb obtained as in (6) 
from the orbitals generated through the orthogonaliza- 

tion of the unperturbed orbitals ! {^fclfeg^j- 

This definition is equivalent to that proposed by Cybul- 
ski and Seversen^^. 

For the proof and the detailed discussion of DFT based 
PB method the Reader is referred elsewhere. 



B. Dispersion- free approximation 

The Eq. (14) shows how the DFT interaction energy 
can be evaluated without refering to supermolecule con- 
cept, using exclusively appropriately perturbed KS equa- 
tions solved for monomers. The decomposition of the 
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DFT interaction energy allows to modify it in such a 
way that we can eliminate the correlation effects between 
two subsystems, so that we obtain the dispersionless in- 
teraction energy between two systems. To this end, we 
describe the interacting monomers with the full xc po- 
tentials while allowing only the exact exchange potential 
operate between them. Mathematically, this involves the 
replacement of the AE^c term (23) in (14) with the ex- 
change interaction which in terms of one-electron reduced 
density matrices (1-DM) for the closed-shell system reads 
(cf. Ref. 26) 

-Ecxch[PA; Pb] = 

= -U [ PA(ri;r,)pB(r2;r0 ^3^^^3^^^ ^^8) 

The 1-DM resulting from a single Slater determinant is 

PA(ri;r2) = 2^a,(ri)a*(r2) (29) 
ieA 

with the density (6) simply being the diagonal part of 
1-DM, 



PA(r) = PA(r;r). 



(30) 



Thus, the dimer energy bifunctional incorporating the 
exact exchange takes the form 



S'abIpa; Pb] = Ea[pa] + Eb[pb]- 
+ £'cist[PA; Pb] + Ecxch[PA', Pb]- 



(31) 



The idea of reducing the intermolecular potential to 
the exchange-only part has been successfully used in 
the model helium dimer calculations by Hefielmann and 
Jansen^'' and Allen and Tozer^^. Now we perform the 
search of the extremals of the bifunctional (31) with re- 
spect to Pa and pe under the constraint of the mutual 
orthogonality between monomers' occupied orbitals. The 
orthogonality constraint ensures that the density addi- 
tivity condition (11) is maintained and the intersystem 
Pauli exclusion principle is fulfilled. To this end, we 
perform the variational optimization in two steps using 
the PB method (see Refs. 20 and 21): first, the bifunc- 
tional extremal search is performed without the imposi- 
tion of the intermonomer orthogonality constraint, and 
secondly, the penalty operator is added in the result- 
ing iterative scheme. The minimization of the bifunc- 
tional (31) leads to a system of coupled equations for 
optimum orbitals: 

(/KS(r) + S^i«t(r) + = eA,M^) 



where the action of the exchange operator on an arbitrary 
one-electron function x reads 



ir\r)x{r) = -I f ^^i^x(r')dV. (33) 



In the second step of the PB procedure, the iterative 
process of solving Eqs. (32) with the aid of the penalty 
operator is formulated in full analogy to (24). The inter- 
action energy at the nth iteration is obtained upon the 
insertion of the densities calculated with orthogonalized 
orbitals resulting from (26) into (31) and subtracting the 
unperturbed monomer energies: 



«PB[n] ^ 
<^i„t = ^AE 



~[7i\ ~[n\ 

Pa ;Pb 



Ea [pI] - Eb [pI] = 



Pa ;Pb 



E. 



cxch 



Pa ;Pb 



(34) 



In the above equation, the A monomer deformation is 



AEa 



Ea 



Pa 



Ea [p1] , 



(35) 



and analogously for monomer B. Interaction energy at 
the zero iteration is then 



HL 



<^AB [pl;p'k]-E[pl]-E[p''B]. (36) 



Henceforth, the energy (36) will be referred to as the 
dispersion-free HL interaction energy. The definition of 
the analog of HL energy at the DFT level (27), which is 
simply the zero iteration obtained with monomer densi- 
ties unperturbed by the interaction, depends on a par- 
ticular functional, as shown by Cybulski and Seversen^^ 
and by Rajchel et al.^^ The dispersion-free HL interaction 
energy (36), in contrast to other definitions, rigorously 
excludes the dispersion interaction. It contains the well 
defined electrostatic and exchange interaction contribu- 
tions, in this case between unperturbed isolated DFT 
monomers. It is closely related to the first order en- 
ergy in the symmetry-based perturbation theory based 
on SAPT(DFT). 

represents the final PB energy calculated 
with (34) using self-consistent orbitals satisfying (32). 
It results from the mutual electric polarization of DFT 
monomers, and, owing to the Pauli blockade procedure, 
it contains exchange contributions. It is related to the 
induction terms of the SAPT formalism with two im- 
portant advantages over the latter: PB sums all electric 
polarization terms to infinity and accounts for accompa- 
nying exchange effects in a consistent manner within the 
DFT formalism. 

In the original HF-based formulation,^'^ the PB proce- 
dure simply restores supermolecular HF interaction en- 
ergy, and obviously neglects any kind of electron cor- 
relation. For the DFT analog, both monomers are de- 
scribed with the full KS operator, but are coupled us- 
ing HF Coulomb and exchange operators [v^^^ and ■5™'^'^ 
in Eq. (32), respectively] built from KS orbitals. Such an 
approach accounts for intramonomcr local electron cor- 
relation leaving out the intermonomer nonlocal contribu- 
tions. The S'^^ represents then "non-dispersion" part of 
the interaction energy that includes the electrostatic, ex- 
change, and induction components. For rare gas dimers 
it is purely repulsive (see section HI). 
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C. Pauli Blockade Plus Dispersion 



Closed-shell atoms and ions 



The total interaction energy, termed PBD for Pauli 
blockade plus dispersion, is obtained by adding to (o^^ 
the dispersion component obtained at the DFT level of 
theory, either from SAPT or by other accurate tech- 
niques: 



J7' 



PBD 



'int 



E, 



disp- 



(37) 



In this work we used the second-order dispersion compo- 
nents from SAPT-DFT.29.30 



FED 



°int 



E 



(2) 



E 



,(2) 

cxch-disp' 



(38) 



For the sake of comparison we also calculated the SAPT 
interaction energies by its most efficient hybrid version, 
SAPT^, i.e. through the second order and corrected with 
the so called " Hartree-Fock delta term",'^^ 5y 



3HF- 



(39) 



It includes all leading electrostatic, exchange, induction, 
and dispersion contributions, all evaluated at the DFT 
level of theory. It should be stressed, however, that the 
i5hf correction comes from HF rather than DFT calcu- 
lations. It provides a rough approximation of higher 
than second-order induction-with-exchange terms which 
are necessary to correct otherwise divergent perturbation 
expansion. More explicitly, the HF supermolccular inter- 
action energy and (5hf satisfy 



E: 



HF 



ind,resp 



^(10) 
cxch 
,(20) 



exch-ind,rcsp ' "HFj 



(40) 



where E, 



(10) 
elst 



and E' 



(10) 

cxch-ind 



are first-order electro- 



static and exchange SAPT-HF energies, and El 



and E. 



(20) 



(20) 

ind.rcsp 



denote second-order induction and 



^exch-ind,icsp 

exchange-indution SAPT-HF contributions calculated 
within the coupled-HF formalism. ■^^''^'^ In this approach, 
all the third-order dispersion-exchange-induction effects 
and dispersion-induction coupling are neglected. 



III. RESULTS AND DISCUSSION 

To test the efficiency of the PBD approach we per- 
formed calculations for two general classes of non- 
covalent interactions: van der Waals complexes of closed- 
shell atoms and ions and H-bonded systems. We applied 
four DFT functionals of the meta-GGA PBE hierarchy: 
PBREV'^** and PWQl^^ (local exchange plus correlation), 
and PBEO^*^ and BSLYP^'^ (a hybrid of local and exact 
exchange plus correlation). All calculations have been 
carried out with the aug-cc-pVTZ basis sets. 



We performed calculations for several diatomic sys- 
tems composed of closed-shell atoms and ions. The Ar2, 
ArNa"^ and ArCl~ have already been shown in Ref. 1, 
here we provide additional examples of HeLi+, He2, and 
Ne2 which are typical closed shells and which provide 
stringent test as their intramonomer dynamic correlation 
is demanding for electronic structure methods. These re- 
sults are compared with SAPT(5 and benchmark values 
which typically originate from high-level supermolccular 
calculations. 

Overall, our potentials are in good agreement with the 
benchmark curves for all systems under consideration, cf. 
Figs. 1-3. As anticipated, the best performance is for the 
PBEO potential. 

One can compare the PBD interaction potentials with 
the potentials obtained from straightforward supermolcc- 
ular calculations using the same functionals corrected 
for basis-set supersposition error (BSSE). For the sake 
of brevity we illustrate this comparison only for Ne2 
(Fig. 4), but qualitatively the results for the other two 
complexes (HeLi+, He2), as well as for the complexes in 
Ref. 1 (Ar2, ArNa+ and ArCl~) are similar. One can 
see in Fig. 4 that the supermolccular interaction energies 
reproduce neither the accurate benchmarks (they are far 
too shallow) nor the " dispersion- free" part as they fea- 
ture a small attraction and shallow unphysical minima in 
the long range. Therefore, they are not appropriate for 
combining with with the pure dispersion term as is done 
in many "DFT-|-dispersion" approaches. ^'^''^^^'^"'^ By way 
of contrast, the " dispersion- free" PB potentials based on 
the same DFT schemes (cf. Fig. 4) are purely repulsive, 
revealing proper asymptotic exponential behavior, and 
thus can be adequately corrected by adding the disper- 
sion contribution. 



PBD: PBEREV 
PBEO 
PW91 
B3LYP 
SAPT5: PBEO 
benchmark 




FIG. 1: Comparison of PBD and SAFTJ interaction energies 
for He2. Benchmark results are taken from Ref. 42. 
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FIG. 2: Comparison of PBD and SAPTiJ interaction energies 
for Ne2. Benchmark results are taken from Ref. 43. 




FIG. 3: Comparison of PBD and SAPTJ interaction energies 
for HeLi"'". Benchmark results are taken from Ref. 44. 
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FIG. 4: Comparison of PB and supermolecular DFT interac- 
tion energies for Ne2. 



B. Hydrogen-bonded and other molecular 
complexes 

A set of hydrogen bonded systems from Boese et al.^^ , 
Jurecka et al."*^ and Halkier et al."*^ supplemented with 
two typical van der Waals molecular complexes, the 
methane and ethylene dimers, have been used as a testing 
set. The selection covers a wide range of qualitatively dif- 
ferent interactions: predominantly tt-tt, H-bonds with a 
large dispersion component (ammonia and HCl dimers), 
and finally strong H-bonds involving shared proton char- 
acterized by a relatively large induction contribution. 
The PBEO functional and aug-cc-pVTZ basis set have 
been used. 

A comparison between perturbation and supermolec- 
ular results requires a consistent treatment of monomer 
geometry in both approaches. In our comparison we per- 
formed PBD and SAPT(5 calculations between monomers 
in their optimal dimer geometries. The databases of su- 
permolecular interaction potentials'*^'^^ include the ef- 
fects of monomer deformations within interaction en- 
ergies. For a suitable comparison with the PBD and 
SAPT5 data we recomputed the monomer deformation 
effects at the CCSD(T) level from monomer geometries 
of Ref. 45 and subtracted them from the supermolecu- 
lar interaction energies. These energies are reported as 
benchmark in Table I. 

The results are shown in Table I. Both the total inter- 
action energies as well as the components are displayed, 
and compared with calculations by the SAPT method, 
and the benchmark values from appropriate references 
(denoted in the rightmost column of Table I). The total 
induction energy is a sum of induction and exchange- 
induction contributions, 

p(2) _ p(2) p(2) 

-^ind-tot — ^ind -^cxch-ind' \^'^) 

PBD agrees favourably with benchmark results and with 
SAPT(5. In certain classes of interactions (see discussion 
below) PBD yields more attractive interaction energies 
than both benchmark and SAPTJ. 

The results in Table I allow one to compare the total 
induction effect in the PB method with that in SAPT5. 
To do this, (fj^f can be compared with the sum: eI^^^_^^^ + 
(5hf- The PB induction is consistently more attractive 
than that from the SAPT(5. The discrepancies appear 
to be the largest for the strong H-bonds with a shared 
proton, such as 0H~-H2 0, as well as ammonia- hydrogen 
halides where the induction effect leads to a significant 
stretching of the proton donor. 

Table II allows us to investigate the need for the 
asymptotic correction (AC). A contribution of AC to PB 
and SAPT energy components, 

SE = • 100%, (42) 

where E^c and E denote energies calculated with and 
without AC, is shown in Table I. For all calculations 
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throughout the Paper the correction scheme of Gruning 
et al."** have been used. It is interesting that (o^^ is much 
less sensitive than E^^'' to the presence of AC. This can 
be explained in view of the fact that both S'^^ and E'-^^ 
involve the mutual cancellation of the electrostatic and 
exchange effects of which only the exchange may have a 
wrong asymptotic behaviour. This attests a robust char- 
acter of (fi^t^. 

Another test has been carried out for the bending po- 
tential of the H2O-HF complex. The water and HF ge- 
ometries and the intersystem distance were taken from 
Ref. 47. The potential V{ip), where if is the angle be- 
tween water plane and the line connecting oxygen and 
HF hydrogen atoms, is the interaction energy scaled so 
that it is zero at the minimum. As can be seen in Fig. -5, 
the PBD method gives a remarkably similar potential to 
CCSD(T), while the SAPT(S underestimates the barrier 
by about 40 %. 




130 140 150 160 170 180 
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FIG. 5: Comparison of interaction energies for H2O-HF dimer 
with respect to the bending angle. 



of systems ranging from noble gas dimers to different, 
weakly and strongly H-bonded dimers. 

The new formalism is based on the concept of interact- 
ing separated monomers. The monomers are described 
within the DFT formalism. Their densities interact un- 
der the constraint of the antisymmetry principle and un- 
der the exact exchange intermolecular potential until self- 
consistency is reached. The resulting interaction energy 
represents the dispersion-free part of the total interac- 
tion energy. The new formalism provides a consistent 
definition of the " non-dispersion" part of the interaction 
energy at the DFT level of theory. The total interac- 
tion energy is obtained by a posteriori adding the DFT 
dispersion contribution from SAPT or other formalisms. 
It should be stressed that the PB treatment does not re- 
quire any extra empirical and adjustable parameters (be- 
sides those that are already used by the DFT description 
of monomers). 

It has been demonstrated that for rare gas dimers, hy- 
drocarbon dimers, and both weak and strong H-bonded 
dimers including ionic interaction, the PB combined with 
monomer description of PBEO provides results that agree 
well with the benchmark values. 

The PB technique may be used for clusters of atoms 
and/or molecules and is not restricted to closed-shell sys- 
tems. It may be combined with different functionals to 
describe the monomers. Similarly, it may be combined 
with a variety of formalisms to calculate the dispersion 
part. 
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TABLE II: Comparison of asymptotic corrections (in %) for PB and SAPT components of interaction energies. 
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